Energy-Water Flows from the IM3 GCAM-USA Scenarios

For the IWPR Q1 Update (Jan 2025)

• Goal: plot a water sankey, energy sankey, energy-water sankey for Q1 update of the EW-Flows project

Energy-Water Sankey

# load the data
water_source_target <- read_csv(paste0("../", data_dir, "water_source_target.csv")) 
allenergy_source_target <- read_csv(paste0("../", data_dir, "allenergy_source_target.csv")) 

allenergy_source_target

unique(water_source_target$source)

 [1] "Arkansas White Red"      "California River"        "Caribbean"               "desalination"            "Great"                  
 [6] "groundwater"             "Lower Colorado River"    "Lower Mississippi River" "Mid Atlantic"            "Missouri River"         
[11] "New England"             "Ohio River"              "Pacific Northwest"       "Rio Grande River"        "runoff"                 
[16] "South Atlantic Gulf"     "Tennessee River"         "Texas Gulf Coast"        "Upper Colorado River"    "Upper Mississippi"      
[21] "water_td_an_W"           "water_td_dom_W"          "water_td_elec_W"         "water_td_ind_W"          "water_td_irr_W"         
[26] "water_td_pri_W"         

unique(water_source_target$target)

 [1] "desalination"                  "groundwater"                   "runoff"                        "water_td_an_W"                
 [5] "water_td_dom_W"                "water_td_elec_W"               "water_td_ind_W"                "water_td_irr_W"               
 [9] "water_td_pri_W"                "Beef"                          "Dairy"                         "Pork"                         
[13] "Poultry"                       "SheepGoat"                     "municipal water"               "biomass (conv CCS)"           
[17] "biomass (conv)"                "biomass (IGCC CCS)"            "biomass (IGCC)"                "coal (conv pul CCS)"          
[21] "coal (conv pul)"               "coal (IGCC CCS)"               "coal (IGCC)"                   "CSP"                          
[25] "CSP (dry_hybrid)"              "gas (CC CCS)"                  "gas (CC)"                      "gas (CT)"                     
[29] "gas (steam)"                   "Gen_II_LWR"                    "Gen_III"                       "geothermal"                   
[33] "hydro"                         "PV"                            "refined liquids (CC)"          "refined liquids (CT)"         
[37] "refined liquids (steam)"       "industry"                      "biomass"                       "Corn"                         
[41] "FiberCrop"                     "FodderHerb"                    "MiscCrop"                      "OilCrop"                      
[45] "OtherGrain"                    "PalmFruit"                     "Rice"                          "RootTuber"                    
[49] "SugarCrop"                     "Wheat"                         "nuclearFuelGenII"              "nuclearFuelGenIII"            
[53] "regional coal"                 "regional natural gas"          "regional oil"                  "unconventional oil production"

# write_csv(water_source_target %>% select(source) %>% distinct(), paste0("../", data_dir, "mappings/source_mapping_w.csv"))
# write_csv(water_source_target %>% select(target) %>% distinct(), paste0("../", data_dir, "mappings/target_mapping_w.csv"))

source_mapping_w <- read_csv(paste0("../", data_dir, "mappings/source_mapping_w.csv"))
target_mapping_w <- read_csv(paste0("../", data_dir, "mappings/target_mapping_w.csv"))

water_source_target_mapped <- water_source_target %>% 
  left_join(source_mapping_w, by = c("source")) %>% filter(!is.na(smap)) %>% select(-source, source = smap) %>% 
  left_join(target_mapping_w, by = c("target")) %>% filter(!is.na(tmap)) %>% select(-target, target = tmap) %>%
  group_by(scenario, source, target, year, units) %>% summarise(value = sum(value))

plot_sankey(water_source_target_mapped)

energy_water_source_target <- allenergy_source_target %>% 
  rbind(water_source_target_mapped %>% mutate(value = value * 1e-1)) %>% 
  complete(scenario, year, nesting(source, target), fill = list(value = 0)) %>% 
  select(scenario, source, target, year, value, units)
  

plot_sankey(energy_water_source_target, "Energy-Water Flows in the US")

write_csv(energy_water_source_target, paste0("../", data_dir, "energy_water_source_target.csv"))